Diffusion of Inhomogeneous Vortex Tangle and Decay of Superfluid Turbulence 
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The theory describing the evolution of inhomogeneous vortex tangle at zero temperature is de- 
veloped on the bases of kinetics of merging and splitting vortex loops. Vortex loops composing the 
vortex tangle can move as a whole with some drift velocity depending on their structure and their 
length. The flux of length, energy, momentum etc. executed by the moving vortex loops takes a 
place. Situation here is exactly the same as in usual classical kinetic theory with the difference 
that the "carriers" of various physical quantities are not the point particles, but extended objects 
(vortex loops), which possess an infinite number of degrees of freedom with very involved dynamics. 
We offer to fulfill investigation basing on supposition that vortex loops have a Brownian structure 
with the only degree of freedom, namely, lengths of loops I. This conception allows us to study 
dynamics of the vortex tangle on the basis of the kinetic equation for the distribution function 
nil, t) of the density of a loop in the space of their lengths. Imposing the coordinate dependence on 
the distribution function n{l,r,t) and modifying the "kinetic" equation with regard to inhomoge- 
neous situation, we are able to investigate various problem on the transport processes in superfluid 
turbulence. In this paper we derive relation for the flux of the vortex line density £{x,t). The 
correspoding evolution of quantity C{x, t) obeys the diffusion type equation as it can be expected 
from dimensional analysis. The according diffusion coefficient is evaluated from calculation of the 
(size dependent) free path of the vortex loops. We use this equation to describe the decay of the 
vortex tangle at very low temperature. We compare that solution with recent experiments on decay 
of the superfluid turbulence. 
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INTRODUCTION AND SCIENTIFIC 
BACKGROUND 

Idea that inhomogeneous vortex tangle evolves in a 
diffusive-like manner appeared quite long ago. Thus, the 
authors of paper [1], who observed the regions of high 
vortex line densities £(r,i) — "plugs" in the channel with 
the counterflowing He II proposed that this phenomenon 
appeared due to diffusion of quantity £(r,t). An attempt 
to describe theoretically these processes was made by [1], 
and [2] , who proposed to introduce the term proportional 
V^C{Y,t) into the Vinen equation. Authors were not able 
to restore the value of the diffusion coefficient from ex- 
periment (these results are reviewed and discussed in [3]). 
In series of works (see [4] and references therein) both 
the diffusion process and the diffusion coefficient we ob- 
tained from the so called Non Equilibrium Thermody- 
namics principles developed by the authors earlier. In 
paper [5] the spatial diffusion of an inhomogeneous vor- 
tex tangle had been studied numerically (see Fig. 1). 
Analyzing their results the authors determined the dif- 
fusion constant to be of order of 0.1k (k is the quantum 
of circulation). Dynamics of the inhomogeneous vortex 
tangle had been studied numerically also in paper [6], 
however, since the authors studied the dilute tangle they 
observed the " ballistic" regime rather then pure diffusion. 
Especial interest to the diffusion processes arises in con- 
text of the problem on decay of the vortex tangle at zero 
(extremely low) temperature. Indeed, the most apparent 
mechanism of dissipation in quantum fluids - the mutual 
friction between the vortices and the normal component 



-disappears. At the same time in a series of experimen- 
tal works [7], [8], [9] the decay of quantum turbulence is 
observed, and the question arises what is the mechanism 
for dissipation at zero temperature. Various approaches 
and ideas such as a cascade like break down of the loops, 
Kelvin waves cascade, acoustic radiation, reconnection 
loss etc. have been discussed in details in recent review 
[10]. The mechanism of the vortex tangle speading with 
the subsequent degeneration is usually ignored. In fact, 
contibution of diffusion had been discussed in the cited 
experimental work, but grounding on the value 0.1k for 
the diffusion constant, the authors concluded that this 
small diffusion coefficient did not lead to correct time 
of decay. In the present paper we develop the theory de- 
scribing the evolution of an inhomogeneous vortex tangle 
on the bases of kinetics of the merging and breaking down 
vortex loops. We showed that evolution of a weakly in- 
homogeneous vortex tangle obeys the diffusion equation 
with the coefficient equal to adout 2.2k, which exceeds 
approximately twentyfold as large the value obtained in 
[5]. We present arguments that the diffusion constant 
would be significantly underestimated in [5] due to espe- 
cial procedure used by the authors. We used the diffu- 
sion equation to describe the decay of the vortex tangle 
at very low temperature. Comparison with the recent ex- 
periments on decay of the superfluid turbulence [7] , [8] , [9] 
is made. 
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FIG. 1: Diffusion of a vortex tangle at t=0 sec(a), t=10.0 
sec(b), t=20.0 sec(c) and t=30.0 sec(d),(Tsubota et, al. [5]. 



THE EVOLUTION OF INHOMOGENEOUS 
VORTEX TANGLE 

Free path 

Vortex loops composing the vortex tangle can move as 
a whole with some drift velocity Vi depending on their 
structure and their length. The flux of the line length, 
energy, momentum etc., executed by the moving vortex 
loops takes place. In the case of inhomogeneous vortex 
tangle the net flux due to gradient of concentration of the 
vortex line density appears. Situation here is exactly the 
same as in classical kinetic theory with the difference that 
the " carriers" are not the point particles but the extended 
objects (vortex loops), which possess an infinite number 
of degrees of freedom with very involved dynamics. In 
addition, while collision (or self-intersection) of elements 
of filaments the so called reconnection of the lines occurs, 
and loops either merge or split, losing their individuality 
and turning into other loops. The full statement of this 
problem requires some analog of the secondary quanti- 
zation method for extended objects, or the string field 
theory, the problem of incredible complexity. Clearly, 
this problem can be hardly resolved in the nearest fu- 
ture. Some approach crucially reducing the number of 
degrees of freedom is required. 

We offer to fulfill investigation basing on supposition 
that vortex loops have the Brownian structure (see [14]). 
In accordance with this approach the average loop can be 
imagined as consisting of many arches with the mean ra- 
dius of curvature equal randomly (but smoothly) con- 
nected to each other. Quantity is important param- 
eter of the approach. It plays a role of the " elementary 
step" in the theory of polymer. It is low cut-off of the 
approach developed, the theory does not describe scales 
smaller then ^o- Statistics of Brownian loops is described 



by the generalized Wiener distribution with the only pa- 
rameter ^0 for all loops. Thus we consider the vortex 
tangle as a collection of vortex loops of various lengths 
I, which is only degree of freedom. Let us introduce dis- 
tribution function n{l,t) of the density of a loop in the 
"space" of their lengths. It is defined as the number of 
loops (per unit volume) with lengths lying between I and 
l + dl. Knowing quantity n{l, t) and statistics of each per- 
sonal loop we are able to evaluate various properties of 
real vortex tangle. The distribution fimction n{l, t) obeys 
the Boltzmann type "kinetic" equation. Study of exact 
solution to this "kinetic" equation allowed us to develop 
a theory of superfluid turbulence, which quantitatively 
describes the main features of this phenomenon [11], [12]. 

This approach turns out to be useful for the study of 
inhomogeneous vortex tangle. In this case we have to im- 
pose the coordinate dependence on the distribution func- 
tion and on parameter ^o, that is to put n{l, r, t), ^o(r, t), 
and to modify the "kinetic" equation with regard to in- 
homogeneous situation. In fact, in this work we restrict 
ourselves to a bit more modest problem, namely we study 
the question of the spatial and temporal evolution of the 
vortex line density C{r,t). The corresponding theory can 
be developed in spirit of the classical kinetic theory with 
the difference that the transport processes are executed 
with the extended objects - vortex loops. Accordingly 
the key questions is to evaluate the drift velocity VJ and 
the free path for the loop of size I. The drift velocity Vi 
is defined via an averaged quadratic velocity of the line 
elements (simple average velocity vanishes due to sym- 
metry) 




The averaging () should be done with use of the Gaus- 
sian model briefly described above [14]. The result of 
according calculations is that Vi = P/y/l^. Quantity 
/? is C„(K;/47r) ln(//ao), where oq is the core radius and 
Cvi is a constant about unity. Velocity Vi can be also 
estimated from the following qualitative consideration. 
Consequently considering the average loop as consisting 
of n = ?/^o arches with the mean radius of curvature 
equal to randomly (but smoothly) connected to each 
other, we take its velocity as the resulting velocity of all 
arches composing the loop. Since the arches randomly 
connected to each other, and have the velocity as for 
rings, Varch = /3/^o (directed along the normal), the re- 
sulting averaged velocity is the "random walking" aver- 
age, = ^^Varch = 

The drift motion is realized until the loops collide with 
other loops with subsequent formation of larger loops. 
Number of collisions Pcoi {dt) per small interval dt can be 
estimated from the "kinetic equation" for the distribution 
fimction n{l) of density of loops in space of their lengths 
I (see [11], [12]). The rate of change of density n{l) due 
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to collisions is 



dn{l,t) 
dt 



= -2 



A{h,iMm2-h' 



■ l)n(l)n,{li)dlidl2. 

(2) 

We omit the processes of the loops breakdown due to 
the self-intersection. The reason is that the migration 
of loops is performed mainly by small loops, the large 
ones undergo reconnections without any essential drift. 
The scattering cross-section A{li, I, I2) describes the rate 
(number of events per unit volume and unit time) of col- 
lision of two loops with lengths I and I2 and forming the 
loop of length li + I = I2. It is evaluated in papers [13] 
A(li, I, I2) — b„iVilil. Here b„i is the numerical factor ap- 
proximately equal to &m ~ 0.2 and Vi = /3/^o is the 
velocity of approaching of vortex line elements (similar 
result for quantity ^ 0.1 0.2 was obtained qualita- 
tively in the context of cosmic strings [16]). Then prob- 
ability Pcoi{dt) for loop to collide with other loops (and 
reconnect) in small interval dt is: 



Pcoiidt) = Adt, (3) 
where quantity A is evaluated with the use of relation (2) 



A = 2 



A{h,l,l2)S{l2 - li - l)n{h)dhdl2. (4) 



We calculate the collision probability A with the use 
of the distribution function n(Z), obtained early n{l) — 
(jl-5/2 ^ggg [11], [12]). Here coefficient C can be ascer- 
tained from the normalization condition that the vortex 
line density C is just C = J n{l)ldl (we used while inte- 
gration the fact that quantity ^0 serves as the low cut-off 
). Simple calculation leads to relation: 



A(/) = 2/36,„£. 



(5) 



In the usual way we conclude that probability P{x) for 
the loop of length / to fly distance x without collision is 



P{x) = 2lCeii-p{~2lbrrXx). 



(6) 



It can be seen from relation (6) that the free path for 
loop of length I is l/2lbm^. 



The flux of length and the diffusion equation 

Knowing the I and ^o-dependent averaged velocity of 
loops and the probability P{x) we can evaluate the spa- 
cial flux of the vortex line density C executed by the 
loops. The procedure is very close to the one made in 
the classical kinetic theory with the difference that the 
carriers have different sizes, this requires additional inte- 
gration over the loop lengths. 




FIG. 2: The net flux of length through the small area element 
placed in a; = and orientated perpendicularly to axis x 



Let us consider the small area element placed at some 
point X and orientated perpendicularly to axis x (See 
Fig. 2). The x-component flux of the line length exe- 
cuted by loops of sizes / placed in 0, tp direction (from 
the left and right sides, correspondingly) and remote 
from the area element at distance R can be written as 
dj^{e,ip,l,R) = I n{e,ip, I, R)Vi cose P[R). The net is 



(dj+((?,(^,/,i?)-dj_(0,^,?,i?)) 



sm 



47r 



dldR. 
(7) 

Taking that spacial density of the loops is the function 
of coordinate x (i.e. n{l,x)) we obtain after simple inte- 
gration 



24 £6„ 



1 



2il'^dx^''^ ildx'^ 



1 d 



(8) 



There are possible the following three variants: i). The 
mean radius of curvature is the fixed quantity (this 
case can be realized while injecting the loops (rings) of 
some fixed size into volume.), ii). Quantity ^0(2;) is in- 
dependent variable, which can appear, when the own 
vortex line (kelvin wave propagation) dynamics prevails 
collisions of loop, determining thereby the structure of 
the vortex tangle. In this case the mean radius of cur- 
vature ^0 should be connected to spectrum of Kelvin 
waves[17] iii). Mean curvature £,i){x) and vortex line din- 
sity C{x) are not independent, the relation between these 
two quantaties appears due to the balancing of fusion and 
splitting of the loops and can be determined from the 
"kinetic" equation [11], it is ^0 « 0.27£-i/2. The fact 
that the interline space is of the order of the mean ra- 
dius of curvature had been firstly discovered numerically 
in [18]. Using this connection we come to conclusion 
that Jx — ~DydC/dx and correspondingly, the spatial- 
temporal evolution of quantity C obeys the diffusion type 
equation 

ac 

'dt 



(9) 



where the diffusion coefficient Dy is equal C^k. Our 
approach is fairly crude to claim for good quantitative 
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description. However, if to adopt the data grounded on 
exact solution to the Boltzmann type "kinetic" equation 
([12]) we conclude that ~ 2.2. As it was discussed 
above, the result that evolution of VLD obeys the diffu- 
sion type equation with the diffusion constant of order of 
K is expected. The actual interest is connected to prefac- 
tor Cd- Early, the only quantitative result was obtained 
numerically by Tsubota et al. , and it was 0.1, in 

22 times smaller. We will discuss the probable reason for 
this large discrepancy later. 



BOUNDARY CONDITIONS 

To construct the boundary conditions we have to state 
the problem specifically. 

1. Smearing of the tangle. Let us consider, first, 
the case when the vortex tangle is placed in some re- 
stricted domain of superfluid helium. Let us consider also 
that vortex line density C is not too high, this allows the 
relatively large loop to be radiated. They move slowly, 
the smaller loops run down larger loop, then collide and 
reconnect with them. So, outside of initial domain the 
well developed tangle is formed. This, secondary, vortex 
tangle smoothly joins the initial tangle inside the domain. 
This implies that in this case no boundary conditions are 
required at all, and evolution of the vortex line density 
obeys equation (9) in infinite space with initial distribu- 
tion £(r,0) inside this domain. 

2. Radiation of loops. The second situation can be 
realized when the radiated vortex loops run away does 
not influence the initial vortex tangle. It can happen, for 
instance if the initial tangle is very dense, so it can radiate 
only very small loops, which propagate rapidly. These 
loops run away without interaction with each other and 
with initial tangle where they are radiated from. Other 
hypothetic variant appears when there is some trap on 
the boundary absorbing vortex loops. Thus the vortex 
loops escaped from the initial domain do not influence 
on the original vortex tangle. In both case the boundary 
conditions can be found assuming that diffusive like flux 
of length near boundary J = — Z)^V£(xf,, t) coincides 
with the flux executed by vortex loops radiated through 
the (right) boundary 3r{xb, t). The latter is evaluated by 
small modernization of (7) 

Jr= [ ln{l, Xb + -Rcos e){vi cos 0)P{R) ^''^''"^ dldR. 
J 47r 

(10) 

Contribution due to remote loops (term Rcos 6 in ar- 
gument for n(l, Xb + R cos 9)) is evaluated as early giving 
obvious result J = —{l/2)Dy\/£{xb,t). To evaluate con- 
tribution from the first one Jgvap we use again solution 
n{l) = Cl~^^'^ and relation (6). Simple calculation lead 
to result that flux 3 evap of length due to evaporation (ra- 



diation) is 

3evap = Cradf^^^ , (H) 

with Crad ~ 0.47. The sum of and Jevap gives the 
flux through the boundary, which has to be equal to 
—Dv'VC{xb,t) at the boimdary. Then we finally have 
the following boundary condition 

CradCi + ^D^VjC{Xb, t) = 0. (12) 

3. Solid walls. In case of solid wall which corre- 
sponds to some experiments the situation is more in- 
volved. Vortices can annihilate on the solid wall, they can 
undergo pinning and depinning radiating vortices back to 
the bulk of helium. Surely this requires a special treat- 
ment which goes beyond the scope of the work. One 
the possible ways is to consider the solid wall as a "par- 
tial" trap, which catches the loops and re-emits a part of 
them back into the volume. Formally it can be written as 
condition (12) with additional term Jback describing the 
back fiux. Without detailed analysis it can be supposed 
that the back flux is proportional to the vortex line den- 
sity on boundary Jback = —Cback^{xb,t) with coefficient 
Cfeocfc depending on dynamics of line on the wall (jumps 
between pinning sites, Kelvin waves dynamics near the 
wall etc.). Thus the boundary condition can be written 
in form. 

CradCi + \D^VC{xb, t) - CbackC = 0. (13) 

DECAY OF THE VORTEX TANGLE 

Let us discuss the question how the approach devel- 
oped can be applied to the problem of decay of the vortex 
tangle at zero temperature. As it was discussed in the 
Introduction, the contribution of diffusion (or radiation 
of loops) is usually ignored, mainly due to smallness of 
the diffusion constant offered in paper by Tsubota et all 
[5] . Let us examine this work more thoroughly. In paper 
[5] there was studied one-dimensional evolution (spacial 
spreading) of the vortex tangle concentrated initially in 
some domain of space and having nonuniform distribu- 
tion there as it is depicted in the first picture in Fig. 
3. 

In the rest pictures distribution of the VLD at differ- 
ent moments of time is shown. To describe this evolution 
of the VLD it had been supposed that quantity £{x,t) 
obeys equation (9) with additional term — X2('«/27r)£^ in 
the right hand side. In turn this term was introduced to 
describe the decay of the vortex tangle in previous nu- 
merical simulation made by Tsubota et al. [15]. Because 
of this additional term contribution of diffusion to the 
whole decay would be significantly underestimated. We 
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FIG. 5; Comparison of experimental data on the decay of 
superfluid turbulence obtained in [8] (upper figure) and our 
theoretical result (lower). We calculated temporal evolution 
of the averaged vortex line density (for initial condition £ — 
10* 1/cm^) due to diffusion process described in the present 
paper. The initial domain of high vortex line density was 
created in the volume "^He-B, so its diffusion like behavior 
should satisfy the first type boundary condition. The straight 
line in the lower figure exactly corresponds to line A in the 
upper figure 5 (which was named by authors of [8] as "limiting 
behavior"). 




0.5 10 O,.'* 1 



FIG. 4: Evolution of the vortex line density calculated with 
use of equation (9) without ani auxiliary term, the diffusion 
constant is equal to Cd ~ 2.2 * 10~^ cm? / s. 

would like to note that there is possible to choose an- 
other reasoning, namely, to consider that decay of the 
vortex tangle in both cited works occurs mainly by the 
diffusion process (with the diffusion coefficient calculated 
in the present work). We calculated spatial -temporal 
evolution of vortex tangle (under condition of numerical 
experiment by Tsubota et al.) with the use of equation 
(9) Results depicted in Fig. 4. enables us to conclude 
that the approach developed describes satisfactory the 
evolution of vortex tangle without any additional suppo- 
sition. 

Thus, we demonstrated that evolution of the vortex 
line density in numerical experiment [15] can be described 
in terms of pure diffusion of the vortex tangle with the 
diffusion coefficient Cd ~ 2.2 * 10"'^ cm? / s. 

Let us now discuss two recent experiments on decay 



of the voretx tangle at very small temperature [8] and 
[9]. The authors reported attenuation of vortex line 
density in superfluid turbulent helium, ('^He-B in paper 
[8] and '*He in paper [9]). They attribute the decay of 
the vortex tangle to the classical turbulence mechanisms. 
Without discussion of this variant we would like just to 
estimate the contribution into attenuation of the vortex 
line density due to the pure diffusion mechanism. 

In the upper picture of Fig. 5 we displayed the Fig 
2. of work [8] showing results of measurements on the 
temporal behavior of the average vortex line density C{t) 
(solid curves, see for details [8]). In the lower picture of 
Fig. 5 we depicted the evolution the same quantity (for 
initial condition C = 10^ l/cvn?) due to diffusion pro- 
cess described in the present paper. The initial domain 
of high vortex line density was created in the volume 
■^He-B, so its diffusion behavior should satisfy the first 
type boundary condition. The straight line in the lower 
picture exactly corresponds to line A in the upper one 
(which was named by authors of [8] as "limiting behav- 
ior"). It is easy to see that there is perfect agreement 
between experimental data and theoretical predictions. 

In contrast to work [8] in the paper [9] the decay of 
vortex tangle in He-II was observed in the closed cube 
with solid walls. In the upper picture of Fig. 6 there is 
depicted the temporal behavior of the average vortex line 
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FIG. 6: Comparison of experimental data on the decay of 
superfluid turbulence obtained in [9] (upper figure) and our 
tlieoretical result (lower). We calculated temporal evolu- 
tion of the averaged vortex line density (for initial condition 
C — 10* l/cm?) due to diffusion process described in the 
present paper. The superfluid turbulence in [9] was created 
in the cubic container with solid boundaries. Therefore we 
had chosen the third type boundary condition with the flt- 
ting parameter Ciack ~ 0.9. 



density Cav{t)- 

We calculated the same dependence on the base of dif- 
fusion equation (9) , result is shown in the lower picture of 
Fig. 6. As discussed above we can ( in the frame of the 
present paper) determine the boundary condition only 
up to coefficient of re-emission Cback ■ Considering it as 
a fitting parameter we have chosen the value of Cback ~ 
0.9. It can be seen that the decay of the vortex tan- 
gle due to diffusion reproduces some feature observed in 
experiment. 

In particular, there is some plateau anticipatory decay 
of the tangle. Full decay of the tangle occurs in time of 
the order of few thousands seconds as it was observed in 
experiments. The slope of curve in interval of most inten- 
sive decrease shows the dependence close to ^ f^^"^ (The 
straight line, its position coinsides with position of the 
straight line in upper picture). Thus, there is again very 
good agreement between experimental data and theoret- 
ical predictions. There is one very interesting by-product 
of the consideration exposed above. It is easy to no- 
tice that for small values of vortex line density the bend 
appears on curve >C(t). It corresponds to the fact that 
the diffusive-like flux of the length vanishes (because of 
vanishing the gradient V£(xf,,i)) and the only first and 
third terms in the boundary condition (13) survive. Let 
us recall that these terms correspond to radiation of loops 
from the bulk to the boundary and to the re-emission of 
loops from the wall into bulk of helium. Comparing this 



terms (for the chosen value Chack ~ 0.9) we obtain that 
equilibrium is reached for value of the vortex line density 
C, of order 50 100 1/cm^. This value can be consid- 
ered as a "background" value of pre-existing vortices in 
helium 

CONCLUSION 

In summary, the theory describing evolution of inho- 
mogeneous vortex tangle at zero temperature was devel- 
oped on the bases of kinetics of merging and splitting 
vortex loops. Using the Gaussian model for vortex loops 
we calculated the (size dependent) free path and mean 
quadratic velocity of vortex loops. With the use of these 
quantities we calculated the flux of the vortex line den- 
sity L(x^€) in inhomogeneous vortex tangle and demon- 
strated that under certain circumstances it satisfies to 
the diffusion like equation with the coefficient equal ap- 
proximately to 2.2k. We used this equation to describe 
the decay of the vortex tangle at very low temperature. 
We compare solution with the recent experiments on de- 
cay of the superfluid turbulence. The good agreement 
with the experimental data allowed us to conclude that 
the diffusion processes give the significant contribution in 
the free decay of the vortex tangle at absence of normal 
component. 
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